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Millstone  Quarry,  Figure  1,  situated  on  Millstone  Point  in  Waterford,  Conn- 
ecticut, was  operated  for  the  extraction  of  granite  construction  materials  from  the 
late  1700's  to  the  early  1950's.  In  1943  the  lower  reaches  of  the  quarry  were  aban- 
doned and  subsequently  began  to  fill  with  brackish  water,  the  salt  being  supplied 
by  seepage  along  the  southern  and  western  walls.  From  1943  to  the  final  operations, 
rock  was  quarried  above  the  rising  water  line. 

The  manner  in  which  the  quarry  developed  its  present  hydrological  conditions 
is  uncertain.  However,  it  is  certain  that  Millstone  Quarry  has  characteristics  not 
evident  in  other  closed  bodies  of  water  within  the  Southeastern  Connecticut  region  (1). 
This  is  one  of  the  reasons  for  establishing  at  Millstone  Quarry  the  USL  Millstone 
Quarry  Research  Facility.  The  facility,  which  performs  underwater  acoustics  test- 
ing, is  housed  on  a barge  roughly  centered  over  the  deepest  portion  of  the  quarry. 

In  order  that  the  sound  velocity  structure  might  be  calculated,  the  vertical  salinity 
profile  wasi determined  and  assumed  to  be  invarient  in  time.  Also  a vertical  string 
of  thermocouples  was  installed  to  obtain  the  time  variations  of  the  vertical  thermal 
profile.  The  temperature  data  discussed  in  this  report  were  obtained  for  just  such 
sound  velocity  calculations  and  were  made  available  through  J.  G.  Navin,  Acting 
Branch  Head,  Millstone  Quarry  and  H.  J.  Wilms,  Acting  Division  Head,  Acoustics 
Measurements  Division,  USL. 

C.  F.  Gerber  (2)  (3)  brought  the  unique  hydrological  conditions  of  Millstone 
Quarry  to  the  attention  of  the  Oceanography  Branch  of  the  Ocean  Sciences  Divi- 
sion, USL,  and  subsequent  investigation  have  resulted  in  detailed  summaries  of  the 
unique  scientific  value  of  Millstone  Quarry  (4)  (5). 
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This  report  is  a first  attempt  at  rationalizing  the  time  variations  of  the  vertical 
thermal  profile  with  existing  theories  on  thermal  mixing  in  the  oceans.  Application 
of  present  theories  on  thermal  mixing  require  evaluating  various  physical  constants. 
This  report  deals  with  the  time  and  space  variations  of  one  of  those  physical  con- 
stants - the  vertical  kinematic  thermal  eddy  diffusitivity,  Kz.  r 


STATEMENT  OF  INVESTIGATION 


When  one  cares  to  make  an  investigation  on  the  inter-relationship  between 
several  inputs  to  the  same  system,  the  obvious  application  of  scientific  principle 
is  to  observe  the  system  output  for  all  possible  controlled  permutations  of  the  several 
inputs.  However,  for  investigations  in  which  environmental  weather  data  are  sys- 
tem inputs,  one  can  not  obtain  in  a controlled  manner  all  possible  permutations  of 
the  several  inputs.  Therefore,  either  exact  measurements  must  be  made  on  each 
of  the  several  environmental  inputs  at  a rate  which  is  greater  than  that  observed 
for  changes  in  the  system  output,  or  assumptions  must  be  made  at  a level  capable 
of  justifying  any  conclusion?  concerning  the  system  output. 


It  is  quite  obvious  that  the  previous  statement  leaves  two  approaches  open 
for  planning  investigations  for  the  generation  of  new  data;  yet  for  any  treatment 
of  historical  data  one  is  forced  to  maleassumptions  about  any  possible  system  inputs 
which  were  not  recorded. 


The  geometry  of  Millstone  Quarry  is  presented  in  Diagram  1 and  since  the  ther- 
mocouples occupy  fixed  points  in  that  geometrical  space,  one  must  assume  that  the 
data  are  adequate  to  completely  represent  the  horizontal  and  vertical  thermal  field. 
Also,  since  temperature  data  were  not  obtained  in  a continuous  manner,  one  must 
assume  that  at  no  time  between  any  two  data  samples  did  the  temperature  have  a 
value  which  is  not  bracketed  by  the  temperature  corresponding  to  the  two  data 
samples.  The  thermocouples  measured  the  temperature  at  a fixed  depth  below  the 
surface,  but  a rise  in  water  level  of  two  feet  occurred  during  the  test  period  and 
caused  a change  in  the  absolute  measurement  geometry.  Therefore  one  must  also 
assume  that  the  vertical  thermal  changes  were  rapid  compared  to  the  rise  in  water 
level. 


It  should  be  noted  that  while  no  data  are  available  to  validate  the  preceding 
assumptions  for  the  period  July  1965  to  August  1966,  investigations  subsequent  to 
that  time  indicate  that  the  assumptions  are  realistic  for  the  period  August  1966  to  the 
present.  It  is  in  this  context  that  figure  2 presents  the  time  variations  of  the  ver- 
tical thermal  profile  for  Millstone  Quarry.  L _ 3 f 5 = * C 
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Diagram  1 Simplified  Geometry  of  Millstone  Quarry 
THEORY  AND  COMPUTATIONS 

Knowledge  is  obtainable  on  ocean  mixing  by  application  of  the  cons..  \ 
laws  of  mass,  momentum,  and  energy  to  the  data  which  one  observes.  The  p.  >»«nt 
discussion  is  confined  to  the  exchange  of  thermal  energy  across  the  Millstone  Quarry 
air/sea  interface  and  its  subsequent  distribution  in  the  vertical  column.  The  con- 
tinual exchange  of  energy  across  the  boundaries  of  the  quarry  require  the  water 
column  to  continually  re-establish  equilibrium  by  internal  redistribution  of  heat  and 
matter. 


For  any  given  set  of  equilibrium  conditions  (temperature,  salinity,  and  pressure) 
one  may  define  the  media  by  the  equation  of  state  for  sea  water.  Small  perturbations 
of  equilibrium  conditions  will  require  changes  in  the  media  to  re-establish  equi- 
librium. Changes  in  the  media  will  be  such  that  the  entropy  of  the  new  equilibrium 
state  will  be  at  a maximum  and  by  examining  the  entropy  of  all  possible  final  states 
having  the  same  total  energy  and  mass,  one  can  define  the  new  equilibrium  con- 
ditions. The  only  additional  information  required  to  describe  the  changing  media 
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is  the  rate  which  the  system  approaches  equilibrium  for  any  specified  change  in 
conditions  (temperature,  salinity,  or  pressure). 

As  previously  stated,  the  data  for  this  report  were  not  obtained  for  the  spe- 
cific purpose  of  determining  the  rates  of  mixing.  The  data  at  hand  are  inadequate 
to  describe  the  transfer  of  salt  within  the  media  except  to  say  that  it  is  small  and 
will  be  disregarded  as  a possible  contribution  to  the  maximizing  of  entropy  (6)  (7). 
Likewise,  pressure  changes  are  neglected  and  the  equation  of  state  is  reduced  to 
changes  in  density  as  a function  of  temperature.  See  Figure  3 (8).  The  salinity 
profile  presented  in  Figure  4 must  be  consulted  to  determine  the  quiescent  density 
for  any  given  temperature  and  depth. 


Heat  flux  at  the  upper  boundary  of  Millstone  Quarry  alters  the  temperature 
at  that  boundary,  and  then  the  change  in  temperature  spreads  into  the  interior. 
The  temperature  changes  spread  away  from  the  boundary  by  a combination  of  mo- 
lecular diffusion  and  turbulent  diffusion.  Turbulent  diffusion,  like  molecular  dif- 
fusion, is  proportional  to  the  product  of  specific  heat  and  vertical  thermal  gradi- 
ent but  generally  with  a much  larger  proportionality  constant  (9). 


Molecular  diffusion  will  always  be  present,  but  turbulent  diffusion  depends 
on  the  presence  of  turbulance.  Turbulence  may  be  generated  at  the  upper  surface 
by  wind  stress  or  may  be  generated  due  to  the  conversion  of  potential  energy  into 
kinetic  energy  as  is  the  case  when  a negative  density  gradient  is  present.  Turbulance 
may  be  decreased  by  the  buoyancy  effect  and  viscosity  in  regions  where  the  density 
gradient  is  positive.  Diagram  2 represents  the  relationships  of  the  preceding  dis- 
cussion. 
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Diagram  2 Energy  flow  diagram  for  Millstone  Quarry  vertical  mixing 
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The  high  quarry  walls  will  reduce  the  effects  of  wind  on  the  upper  surface  (10) 
and  in  the  absence  of  wind  stress  a positive  heat  flux  will  cause  the  system  to  tend 
toward  vertical  stability.  If  the  heat  flux  is  negative/  the  system  will  tend  toward 
vertical  turbulance. 
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The  thermal  spreading  is  represented  as  ^ ^ - V ^ 


(1) 


For  the  quarry  one  must  assume  that  the  substantial  derivative  is  equivalent  to  the 
local  derivative. 


Where 
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0=  CpT 

Cp=  specific  heat  at  constant  pressure 
T = temperature 

Kz=  vertical  exchange  coefficient  for  heat. 


(2) 


If  one  assumes  that  in  time  the  specific  heat  at  any  depth  is  constant  and 
that  Kz  and  C p are  constant  over  the  depth  interval  between  two  measurements, 
one  may  express  the  kinematic  eddy  diffusitivity  as 

JTr 

J JL* 

By  entension  to  Proudman  (11),  one  may  consider  the  determination  of  K»  from  the 
time  variations  of  isotherms  at  a given  station.  See  Diagram  3.  Draw  lines 
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a measurement  was  obtained.  Let  denote  the  increment  of  temperature  sep- 
arating isotherms,  the  increment  of  average  temperature  from  Z to  (Z-  z) 
and  denote  the  increment  of  average  temperature  from  Z to  (Z+  z). 

Let  denote  the  time  separation  between  two  isotherms  at  the  same  depth  and 
let  all  derivatives  be  evaluated  at  P. 


Therefore,  in  time  at  any  given  depth  Z, 

f r ir  j> 

Also  from  the  curve  for  the  average  temperature  as  a function  of  depth 


(4) 
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P!  PP"  is  considered  to  represent  the  average  thermal  variations  within  the 
neighborhood  bordered  by  the  two  isotherms  separated  by  J" jnd  the  two  depth 
limits  separated  by  2JL  The  evaluation  of  equation  6 is  graphically  determined 
as  per  diagram  3. 
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which  is  independent  of  the  temperature  difference  between  the  two  isotherms. 

Kz  as  determined  from  equation  7 are  presented  in  Figure  5. 

RESULTS  AND  CONCLUSIONS 

A comment  on  dimensional  units  is  in  order.  The  isotherms  in  Figure  2 are 
presented  in  the  English  system  and  are  the  units  in  which  the  measurements  were 
taken.  It  has  been  previously  stated  that  the  data  were  not  taken  specifically 
for  this  study.  The  depth  units  on  Figure  4 and  5 were  reported  in  feet  to  allow 
ease  of  comparison  with  Figure  2.  This  is  contrary  to  scientific  convention, 
however,  the  necessary  conversions  were  made  to  allow  Kz  to  be  reported  in  the 
conventional  units  of  cm  2/sec.  The  combinations  of  mixed  units  presented  in 
this  report  in  no  way  detract  from  the  subject  at  hand.  There  was  however  a 
serious  misleading  artifact  introduced  into  Figure  5.  The  decimal  positions  are 
not  indicative  of  the  number  of  significant  figures.  The  values  are  only  accurate 
to  the  nearest  one-tenth  cm  2/sec. 

The  Kz  in  figure  5 are  heavily  averaged  in  time,  first  because  the  data  were 
obtained  at  most  once  a day  which  immediately  places  a lower  limit  on 
Furthermore,  data  were  available  for  only  198  days  of  the  385  days  interval 
covered  in  this  report.  Additional  time  averaging  was  necessary  to  eliminate 
scatter  resulting  from  the  evaluation  of  the  individual  thermal  profiles. 

The  accuracy  of  the  individual  thermal  profiles  is  ±0.5*  and  the  precision 
of  any  point  on  a given  profile  is  ±0.  1°  F. 

For  those  portions  of  Figure  2 where  the  isotherms  approach  vertical  and  Kz 
is  very  large,  the  maximum  Kz  that  can  be  reliably  determined  from  the  figure  is 
1.  2 cm  /sec  for  the  measurements  at  three  foot  separation  and  3.  3 cm  Vsec  for  the 
five  foot  intervals.  For  those  near  vertical  isotherm  portions  of  Figure  2,  one  may 
return  to  the  individual  thermal  profiles  and  determine  Kz  provided  the  tem- 
perature differentials  exceed  the  precision  and  accuracy  of  the  measurement  system. 
At  the  system  limit,  the  Kz  would  be  five  times  that  determined  from  Figure  2. 

For  those  portions  of  Figure  2 where  the  isotherms  approach  horizontal  and  for 
the  region  below  the  56F  isotherm  the  Kz  is  drastically  reduced.  The  minimum  Kz 
that  could  be  detected  from  Figure  2 is  8.6  x 10  ”^cm  2/sec 

There  are  many  combinations  of  comparisons  which  one  might  employ  in 
evaluating  Figure  5.  Only  three  will  be  discussed.  The  time  variations  of  Kz 
in  the  depth  interval  between  33  and  38  feet  serve  to  demonstrate  the  general 
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pattern  of  increasing  during  Fall  and  Winter  and  decreasing  in  Spring  and  Sum- 
mer. This  pattern  is  evident  to  a great  extent  in  the  regions  above  the  halocline. 
The  vertical  variations  of  Kz  are  quite  different  for  the  various  seasons  but  the 
important  observation  is  that  starting  very  rapidly  in  the  Fall  and  continuing 
through  the  Winter,  large  portions  of  the  vertical  water  column  obtain  the  same 
degree  of  turbulance.  This  indicates  that  essentially  the  entire hyperlimnion 
will  become  thoroughly  mixed.  This  is  the  same  conclusion  one  might  have 
obtained  from  the  salinity  profile  of  Figure  4. 

Another  important  feature  of  Figure  5 is  that  during  the  Spring  and  Summer 
the  lower  region  of  the  hyperlimnion  lag  behind  the  upper  regions  in  the  decrease 
of  turbulance.  This  is  because  warming  of  the  surface  layers  which  decreases  the 
turbulance  has  thereby  decreased  the  rate  at  which  heat  can  be  transported  to 
an  underlying  water  mass. 

The  third  point  to  be  discussed  is  the  almost  complete  isolation  of  the  hyp- 
olimnion  from  the  thermal  variations  produced  by  heat  flux  on  the  quarry's  upper 
surface.  This  isolation  is  due  to  the  large  positive  density  gradient  resulting 
from  the  thermal  and  salinity  gradients.  The  buoyancy  effect  is  sufficient  to  damp 
out  turbulance  thereby  preventing  thermal  exchange  between  the  hyperlimnion 
and  hypolimnion.  This  depression  of  turbulance  is  sufficient  to  allow  the  hypo- 
limnion  to  be  anoxic.  (7) 

In  summary  the  energy  relationship  of  Diagram  2 is  substantiated  by  the 
temporal  and  spatial  variations  of  turbulance  as  demonstrated  in  Figure  5. 

Despite  the  limitations  imposed  on  this  report  by  the  nature  of  the  data,  it  is 
clearly  evident  that  within  the  limits  of  Millstone  Quarry  the  vertical  kinematic 
eddy  diffusitivity  varies  over  three  orders  of  magnitude.  This  range  is  in  excess 
of  that  normally  observed  in  the  oceans  and  on  the  lower  end  is  only  six  times  that 
attributed  to  molecular  diffusion.  (12) 
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Fig.  3 Sigma  T for  sea  water  as  a function  of  temperature  and 

salinity. density  maximum; freezing  point. 


Fij.  4 Salinity  Profile  Millstone  Quarry 


(JSL  Tech  Memo  No.  2213-d-67 


J.  S.  Navy  Jnderwater  Sound  Laboratory 
NP24  - N29415  - 12  - 66 


Official  Photograph 


U.  S.  Navy  Underwater  Sound  Laboratory  Official  Photograph 

NP24  - N29652  - 2 - 67 


